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Abstract. Fractional Ornstein-Uhlenbeck process of the second kind 
>— -^ ' (fOU2) is solution of the Langevin equation dXf = —OXt At + dY^ , 6 > 

with driving noise Y/ :— L e~°di3aj,; at — HeTT where B is a 
fractional Brownian motion with Hurst parameter H £ (0, 1). In this 
article, we consider the case H > ^. Then using ergodicity of fOU2, we 
construct consistent estimators of drift parameter d based on discrete 
O^ ' observations in both cases: (i) the Hurst parameter H is known and (ii) 

the Hurst parameter H is unknown. Moreover, using Malliavin calculus 

. - technique, we prove central limit theorems for our estimators which is 

Ph ' valid for the whole range ff G (|, 1). 
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►^ ■ 1. Introduction 

/^ \ Assume B = {-Bt}f>o is a fractional Brownian motion with Hurst parameter 

T-|- ' i7 € (0, 1), i.e a continuous, centered Gaussian process with covariance 

CN ■ function 

^': RH{s,t) = -{s^^ + t^^ -\t-s\''"], s,t>0. 

O ■ 2 

Cn ' Consider the following Langevin equation 

(LI) dXt = -eXt dt + dNt, > 0, 

with some driving Gaussian noise N. When driving noise N = B is fractional 
Brownian motion, solution of the Langevin equation (jl.ip is called fractional 
C^ ■ Ornstein-Uhlenbeck process of the first kind, in short (fOUi). Fractional 

Ornstein-Uhlenbeck process of the second kind is solution of the Langevin 
equation (jl.ip . when driving noise Nt = Y^ = J^ e~^ dBa/, at = HeH . 
The terminologies "of the first kind" and "of the second kind" are taken from 
Kaarakka and Salminen [S] . The motivation behind fOU2 process is that it is 
related to Lamperti transformation of fractional Brownian motion. We give 
more information on fractional Ornstein-Uhlenbeck processes in subsection 

Statistical estimation of drift parameter 6 based on data recorded from con- 
tinuous (discrete) trajectory oi X is an interesting problem in realm of math- 
ematical statistics. In the case of fractional Ornstein-Uhlenbeck process of 
the first kind, the two popular estimators, maximum likelihood (MLE) and 
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2 AZMOODEH AND VIITASAARI 

least squares (LSE) estimators based on continuous observations of the pro- 
cess are considered in Kleptsyna and Breton [TO] and Hu and Nualart [7J 
respectively. In this case, it turns out that MLE and LSE provide strong 
consistent estimators. Moreover their asymptotic normality is proved in [5] 
for MLE and iJ > i and in [7j for LSE and H G [i, |). In the case of 
fractional Ornstein-Uhlenbeck process of the second kind, Azmoodeh and 
Morlanes [2] showed that LSE is a consistent estimator with continuous ob- 
servations. Moreover a central limit theorem for LSE is obtained for H > i^- 

However, from practical point of view, it is natural to assume we have a 
data collected from process X observed at discrete times. In discrete and 
fractional Ornstein-Uhlenbeck process of the first kind case, estimation of 
drift parameter 9 with discretization procedure of integral transform con- 
sidered by Xiao et. al. in [12] assuming the Hurst parameter H is known. 
In the same setup, Brouste and lacus ^4] give an estimation procedure that 
could estimate both drift parameter and Hurst parameter H with discrete 
observations. 

In this article, we take the advantage of ergodicity of fOU2, when > 0, 
to construct consistent estimator of drift parameter 6 based on observations 
of the process in discrete times. Assume that we observe the process in 
discrete times 0, Atv, 2 A at, • • • , A^Aat and Tjq = NAn denotes the length of 
the "observation window". Our aim is to show that: 

(i) when H is known one can construct a strong consistent estimator 6 
with asymptotic normality property under the mesh conditions 

Tn -^ oo, and iVA^ -^ 

with arbitrary mesh A^v such that Ajv — )■ as A^ tends to infinity. 

(ii) w^hen H is unknown one can construct another strong consistent esti- 
mator 9 with asymptotic normality property under restricted mesh condition 

A^ = Ar-, with a G (i, ^^^^ A 1). 

The paper is organized as follows. In section 2, we give auxiliary facts on 
Malliavin calculus and fractional Ornstein-Uhlenbeck processes. Section [3] 
is devoted to estimation of drift parameter when H is known. In section 
m we give a short explanation of estimation of Hurst parameter H based 
on discrete observations. Section [5] deals with estimation of drift parameter 
when H is unknown. We also collect all technical computations in section 
appendix A. 



2. Auxiliary facts 

2.1. A brief review on Malliavin calculus. The use of Malliavin calculus 
to obtain central limit theorem for a sequence of multiple Wiener integrals is 
well established. In this subsection, we briefly introduce some basic facts on 
Malliavin calculus with respect to Gaussian processes. For more details, we 
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refer to [T], [15] and [H]. Let W he a Brownian motion and G = {Gt}te[o,T] 
is a continuous centered Gaussian process of the form 



Gt= I K{t, s)dWs 
Jo 



where the kernel K, < s < t < T satisfying sup(g[o yj /g K{t, s)^ds < oo. 
A typical example of this type of Gaussian processes is fractional Brownian 
motion B with Hurst parameter H G (0,1). It is known that in the case 
H > 2, the kernel takes the form 

1 /"* 3 1 

KH{t,s)=CHS"^ j {u-s)"-2u"-2du. 

Moreover, we have the following inverse relation 

(2.1) Wt = B{iK*H)-\l[o,t]))^ 
and the operator K'^ is defined as 

{K*j,^){s) = l\{t)^{t,s)dt. 

Consider the set £ of all step functions on [0,T]. The Hilbert space Ti 
associated to process G is the closure of £ with respect to inner product 

(l[0,t]'l[0,s])w = RG{t,s), 

where Rcitjs) denotes the covariance function of G. Then, the mapping 
l[o,t] ^ Gt can be extended to an isometry between Hilbert space 11 and 
Gaussian space T-Li associated with Gaussian process G. Consider the space 
S of all smooth random variables of the form 

(2.2) F = f{G{^i),--- ,G{^n)), fi,--- ,Vn€^n, 

where / G C^(IR"). For any smooth random variable F of the form 12.21 we 
define its Malliavin derivative D^ > = D as an element of L'^{Q;'H) by 

71 

DF = Y,^^f{G{ipl),■■■ ,G{ipn))^i. 
i=l 

In particular, DGt = l[o,t]- We denote by Hq = B^'^ the Hilbert space of 
all square integrable Malliavin derivative random variables as the closure of 
the set S of smooth random variables with respect to norm 

\\F\\l2 =JE\F\'^ + JE{\\DF\\'^). 

The divergence operator denoted hy 5g = S is the adjoint operator of Malli- 
avin derivative operator D. For any "H-valued random variable u G L^{Q,;'H) 
belongs to the domain DomJ of the divergence operator 6, then the random 
variable 5{u) is defined by duality relationship 

IE(F5(u)) = IE(DF, u)n, VF € D^'^ 

Also, an element u G L'^{i};T-l) belongs to the domain Dom6 if and only if 

\JE{DF,u)n\<Cu\\F\\L2, 
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for any F G ID^'^, where c„ is just a constant depending on u. The divergence 
operator 6 is also called Skorokhod integral. It is known that in the setup 
of Brownian motion it coincides with Ito integral for adapted integrands. 
Consider the linear operator K* from £ to L^[0, T] defined by 

{K*ip){s) = v{s)K{T, s)+ [ [^{t) - ^{s)] K{dt, s). 

J s 

The Hilbert space Ti generated by covariance function of the Gaussian 
process G can be represented as T-L = {K*)~^{L'^[0,T]) and D^j {71) = 
{K*)-'^{I])]^{L'^[0,T])). For any n > 1, let ^ be the nth Wiener chaos 
of G, i.e. the closed linear subspace of L^(0) generated by the random 
variables {Hn {G{ip)) , if G 7i, \\ip\\H = 1}, and Hn is the nth Hermite poly- 
nomial. It is well known that the mapping /^((/?®") = nlHn {G{(p)) provides 
a linear isometry between the symmetric tensor product ?^®" and subspace 

The next proposition is taken from [16j and provides a central limit theorem 
for a sequence of multiple Wiener integrals. Let AA(0, o^^ denote the Gauss- 
ian distribution with zero mean and variance a^ . The notation — > stands 
for convergence in distribution. 

Proposition 2.1. Let {Fn\n>\ ^e a- sequence of random variables in the q- 
th Wiener chaos, q >2, such that lim„_5.oo ^{Fn) — ^^- Then the following 
statements are equivalent: 

(i): Fn converges in distribution to M{0,a'^) as n tends to infinity. 
(ii): ||L)F„|||^ converges in L^(J7) to qa^ as n tends to infinity. 

2.2. Fractional Ornstein-Uhlenbeck processes. In this subsection, we 
briefly introduce fractional Ornstein-Uhlenbeck processes. The main refer- 
ences are 0, and |2]. We mostly focus on fractional Ornstein-Uhlenbeck 
process of the second kind. Moreover, we provide some new results on frac- 
tional Ornstein-Uhlenbeck process of the second kind. Let B = {Bt}t>o 
be a fractional Brownian motion with Hurst parameter H G (0, 1). To ob- 
tain fractional Ornstein-Uhlenbeck process, consider the following Langevin 
equation 

(2.3) dui^'^°^ = -dU^^'^^Ut + dBt, 

with initial value ^o > The solution of the following SDE can be expressed as 

(2.4) ui'''^'^ = e-''Uo+ f e'^dB, 

Notice that the stochastic integral is understood as a path-wise Riemann- 
Stieltjes integral. Let B denotes a two sided fractional Brownian motion. 
The special selection 
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for the initial value .^o leads to an unique stationary, Gaussian process U^^' 
of the form 

(2.5) [/i^)= r e-^(*-^)d4. 

Definition 2.1. We call the process U^ '^'^' given by ()2.4p a fractional 
Ornstein- Uhlenbeck process of the first kind with initial value (,q . The pro- 
cess U^^' defined in (j2.5p is called stationary fractional Ornstein- Uhlenbeck 
process of the first kind. 

Remark 2.1. It is shown in [5] that the covariance function of the stationary 
process U^^' decays like a power function, so it is ergodic and for H G (2) 1); 
it exhibits long range dependence. 

Now, we define a new stationary, Gaussian process X^°^' by means of Lam- 
perti transformation of fractional Brownian motion B: 

Xi"^ := e-'^'Ba,, t E R, 

where a > and at = "^Ch. We aim to represent the process X^""' as 

solution of a Langevin type SDE. For this reason, consider the process Y°' 

defined via 

pt 

Yi»).= / e-"MB„,, t>0. 
Jo 
The stochastic integral is understood in as path-wise Riemann-Stieltjes in- 
tegral. Using the self-similarity property of fractional Brownian motion one 
can see that the process y'°) satisfies in the following scaling property 

(2.6) {Y};]h>o'= {a~"Y^^}t>o, 

where = stands for equality in law. Using y(") the process X^°'' can be 
viewed as the solution of the following Langevin type SDE 

dX^^ = -aXi''^ dt + dY}''\ 
with random initial value X^"^ = fi^,, = Bn/a ~ ^(0, (f )^^). 

Inspired by the scaling property ()2.6p . we consider the following Langevin 
equation with y^^^ as the driving noise 

(2.7) dXt = -eXt dt + dY^^\ e>o. 
The solution of the SDE (|2.7p is given by 

(2.8) Xt = e-^* (Xo + J^ e'^ dYJ'^^ = e"^* (xo + J^ e^'-^> dB,)j 

with Q = 1 in oj. Notice that the stochastic integral is understood as path- 
wise Riemann-Stieltjes integral. The special selection Xq = /_^ e^^"^)* di?as 
for the initial value Xq leads to an unique stationary, Gaussian process 

(2.9) Ut = e-'' f e^'-^>dBa,. 

J —00 
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Definition 2.2. We call the process X given by (j2.8p a fractional Ornstein- 
Uhlenbeck process of the second kind with initial value Xq. The process U 
defined in (j2.9p is called the stationary fractional Ornstein- Uhlenbeck process 
of the second kind. 

For the rest of the paper we assume Hurst parameter H > i^- 

Proposition 2.2. [9j The covariance function of the stationary process U 
decays exponentially and has short range dependence. More precisely 

c{t) := E{UtUo) = ( exp (^ - mm{6', ~ }t) j , as t -^ oo. 

Let vjj be the variogram of the stationary process U, i.e. 

vuit) := ^E {Ut+s - Usf = c(0) - c(t). 

The foUowing Lemma tehs us the behavior of the variogram function vu 
around zero. We will use this lemma in section HI The B(x,y) stands for 
complete Beta function with parameters x and y. 

Lemma 2.1. The variogram function vjj satisfies 

vu{t) = Ht^^ + r{t), and r{t) = oif^) as t ^ 0+. 

Proof. Using Proposition 3.11 of [9], there exists a constant C{H,9) = 
H{2H - l)H^H{i~e) ^^^^ ^Yiai 

c{t) = C{H,9)e-'' (r Hixyf-'^'^lx - j/l^^-^dxdy] . 

Denote the term inside parentheses by $(t). Then with some direct compu- 
tations, one can see 

an 

$(i) = / r* (xy)(^-i)^(y - xf'-^dxdy 
Jo Jo 

+ r Hxyf-'^'^ix-yf^-'dydx 
Jo Jo 

faQ f"y 

+ / (xy)(^-i)^(y-x)2^-2dxdy 

3 J^y 

29 H 1 /--^ 

I^B{{e - l)H +l,2H-l) + J^(ar - aD [ ' ^^'-'^H^ " zf-'^^z. 



_ % 

OH 

Therefore 



^^^^ ^ (2g p^^'^ BUe - l)H + 1,2H- l)e-^* 



2c' Jo 

= c(0) - {2H - l)H^^ xtx [' z(^-i)^(l - z)2^-2dz + n{t), 

J^ 
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where ri{t) = o{t^^) as t — )■ 0^ and is given by 
{2H - 1)H^^ 



ri{t) = -^ B{{e - l)H + l,2H - 1) ( e"''* - 1 - Ot 



29 . .JO 

Let Bincip,Q,x) and 2Fi{a,b,c;x) stand for incomplete Beta function and 
Gauss Hyper-geometric function respectively (sec [18' for definition). Using 
a relation between these two functions ([18J, page 289), we obtain that 

^' ^(^^-1)^(1 - zf'^-^dz = BW^H -1,(0- l)H, 1 - ^) 
^ at 

TTTT-2H 

= H— r^"^'^'^*K - «o)^^-^ 2F,{{0 + 1)H, 1, 2H; 1 - ^) 
ZH — 1 at 



^-ie+i)t^^2H^i ^ ^^(^)) ^^^((^ ^ 1)^^ 1^ 2F; 1 - ^) 






2i7-l 

where rz{t) = o(t^^^^) as t — )• 0^ and is given by 

rz{t) = e-(^+i)* 2Fi{{e + l)i7, 1, 2H- 1 - ^) r^it) 

at 



+ e-l<^+i)* _ 1 2Fi((0 + l)i/, 1, 2H; 1 - ^) t' 

at 



+ e-^''+^^* ( 2Fiiie + l)H, 1, 2H; 1 " ^) " 1 1 *' 

Hence, it is enough to take r{t) = ri{t) — {2H — 1)H^^ t r^{t) and the claim 
follows immediately. D 



In the general solution (|2.8p . take Xq = 0. Then the corresponding fractional 
Ornstein-Uhlenbeck process of the second kind takes the form 

(2.10) Xt = e-'' [ e^'-^'>'dBa,. 

Jo 

Note that, we have the relation 

(2.11) Ut = Xt + e-''C, i= I e^''^>dBa,. 

The following lemma tells us that the fractional Ornstein-Uhlenbeck process 
of the second kind X given by ()2.10p . can be obtained (in distribution sense) 
with integrating the same function as in fOUi with respect to a Gaussian 
noise violated of fractional Brownian motion with a regular kernel. 

Lemma 2.2. [2] Denote Bt = Bt+n — Bh the shifted fractional Brownian 
motion. Let X he the fractional Ornstein-Uhlenbeck process of the second 
kind given by (|2.10p . Then there exists a regular Volterra-type kernel L such 
that 

(2.12) Xt'= fe-'^'-'^dGs 

Jo 
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where the Gaussian process G is given by 



Gt= f (^KH{t,s) + L{t,s))dWs 



and Brownian motion W is related to B by (|2.ip . 

Remark 2.2. Note that by direct computation and applying Lemma 4-3 of 
[2] , the inner product of Hilbert space T-L generated by covariance function of 
Gaussian process G is given by 

J2H-2 r f^ ^f^^^^i^i^^^^iu+v){±^l) 

Jo Jo 
where an = H(2H — 1). 



(<^,V')^ = Q^i7^^-^ / / ip{u)^{v)e^ 



ea — eH 



2H-2 



dvdu 



The following lemma plays the essential role in the paper. To be more 
precise, we use this lemma to construct our estimator for drift parameter. 

Lemma 2.3. [2] Let X be the fractional Ornstein-Uhlenbeck process of the 
second kind given by (j2.10p . Then as T —)• oo, we have 

^1 X?dt^^i9) 

almost surely and in L^, where 

(Off — '\) ff^^ 

(2.13) ^{9) = ^ -^ B{{9-1)H + 1,2H-1). 

9 

The next lemma studies the regularity of sample paths of fractional Ornstein- 
Uhlenbeck process of the second kind X. The Holder continuity of sample 
paths is the message of the lemma and more information on Holder constants 
is given. See also [9j, Proposition 3.4. 

Lemma 2.4. Let X be the fractional Ornstein-Uhlenbeck process of the sec- 
ond kind given by (|2.1U|) . Then for every interval [S,T] and every e > 0, 
there exist random variables Yi = Yi{H,9), I2 = Y2{H,9,[S,T]), Y^ = 
Y^{H, 9,[S,T]), and Y^ = Y^{H, e, [S, T]) such that for all s,t £ [S, T] 

\Xt -Xs\<iYi + Y2 + Y3) \t-s\+ Y^lt - s\"-' 

almost surely. Moreover, 

(i): Yi < 00 almost surely, 

(ii): Y,iH,9,[S,T])'^YkiH,9,[0,T - S]), k = 2,3, 

(iii): Y^{H, e, [S, T]) "^ Y,{H, e, [0, T - S]). 

Furthermore, all moments of random variables Y2, Y^ and I4 are finite, and 
Y2{H, 9, [0, T]), Y3{H, 9, [0, T]) and Y^iH, e, [0, T]) are increasing in T. 

Proof. Assume that s < t. By change of variables formula, we obtain 



where 



Xt = e Bat ~ 6 Bag — Zt, 



Zt = e-'' fBaJ'-'^-du. 
Jo 
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Therefore 

\Xt - Xs\ < \Ba,\\e-'^' - e-'^'l + e-'\Ba, - BaJ + I^^Jje" 



*-e-1 



+ 



'9t 



Jo Jo 



= h+l2 + h + h. 

For the term Ii, we obtain 

h <0\Bao\\t-s\ 

where ^|-Bao| is ahnost surely finite random variable. For the term J3, we 
get 

/3 < sup e""|BaJ|i- sl- 
ue [5, T] 

Note that Z is a differentiable process. Hence for the term /4, we get 



h< 



' sup \Zu\ + sup e ""l^a^l 

u£[S,T] ue[S,T] 



\t- s\ 



Moreover, by using (j2.1ip . we have 

\Xt\<\Ut\ + \^\. 
As a result, we obtain 

\Zu\<\Uu\ + \^\ + \Ba,\ + \e-''BaJ. 
This implies that 



h< 



' sup \Uu\ + e\^\ + e\Ba,\ + {e + 1) sup e-"|s« 

u£[S,T] ue[S,T] 



t- s . 



Collecting the estimates for /i, I3 and /4, we obtain 



h+h+h< 



+ 



2e\Ba,\ + e\c\ 



t- s 



• sup \Uu\ + {e + 2) sup e-"|Sa 

«6[5,T] uelS,T] 



Put 



Yi = 2e\Ba,\ + e\C\, Y2iH,e,[S,T]):=e sup |t/J 

«e[5,T] 

and finally 

Y3iH,e,[S,T]):={e + 2) sup e-"|i?aj. 

«e[5,T] 

Obviously for the random variable Yi the property (i) fulfills. Notice that Ut 
and e~^Bat are continuous, stationary Gaussian processes. So the property 
(ii) follows. Moreover, all moments of supremum of a continuous Gaussian 
process on a compact interval are finite (see [11] for details on supremum of 
continuous Gaussian process). So, it remains to consider the term I2. By 
Holder continuity of the sample paths of fractional Brownian motion, we 
obtain 

h < e-'C{Lo,H,e,[S,T])\at - asl""-' 
<C{oj,H,e,[S,T])\t-s\"-'. 
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To conclude, we obtain (see [I7| and remark below) that the random variable 

C(w, H, e, [S, T]) has all the moments and C{uj, H, e, [S, T]) ^= C{uj, H, e, [0, T 
S]). Now it is enough to take Y^ = C{uj, H, e, [S, T]). D 

Remark 2.3. ^ITj The exact form of the random variable C{uj,H,e,[0,T]) 
is given by 

(rp rp 2 \ 9^ 

/ / ^^'~^:„ dtds] , 
Jo Jo \t - s\ ^ J 

where Cu.e is a constant. We also have W,C{oj, H,e,[0,T])P < c^^pT"^^ for 
some constant c^^p. 

3. Estimation of drift parameter when H is known 

We start with the fact on invertibility of the function ^ which allows us to 
construct an estimator for drift parameter 6. 

Lemma 3.1. As a function ^ : 1R_|- — t- IR^, ^{&) given by 112. 13\) is bijective, 
and hence invertible. 

Proof. Evidently, for > we have ^{0) > 0. Since Beta function is 
differentiable with respect to its arguments, so the function $ is continuous 
and differentiable in 9. It is also evident that ^{6) converges to infinity as 
6 tends to zero and ^{9) converges to zero as 9 tends to infinity. Hence ^ 
is surjective. Moreover, we have 

(nfj_ 1 \ tf2H 

*'(^) = - a2 B{{9 - 1)H + 1,2H- 1) 

+ ^^^^^^§-^Bii9-l)H + l,2H-l)<0, 

because B{x,y) is decreasing in x. Therefore ^ is strictly decreasing and 
hence injective, which completes the proof. D 

We continue with the following central limit theorem. 

Theorem 3.1. Let X be the fractional Ornstein-Uhlenbeck process of the 
second kind given by ()2.10p . Then an T tends to infinity, we have 

Vf (i £ X^dt - ^{9)^ -^ Af{0, a') 

where the variance a"^ is given by 



(3.1) ^^ i[0,Oo)3 L 

X (l — e~H 



e H — e H 



2H-2- 

dzdxdy. 



The proof relies on two lemmas proved in the appendix where we also show 
that (T^ < oo. The variance o"^ is given as iterated integral over [0,oo)^ and 
the given equation is probably the most compact form. 



Ft = ^liCg), 
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Proof of Theorem \3.1\ For further use, put 

(3.2) 

where the symmetric two variables function g is given by 

g{x,y) 

By Lemma 12.21 we have 

Xt '= If {h{t, •)) , h{t, s) = e-^(*-^)l,<t. 

Using multiphcation formula for multiple Wiener integrals and Fubini's the- 
orem, we infer that 



1 

20 L 



^-e\x-y\ _ ^-e(2T-x-y) 



T 



\\Ht,-)\\ldt + ^li 



{h{t,-)^h{t,-))dt 



^^ ^xfdt + ^iiCg)- 



we get 



(3.3) Vr(i^ x^dt-^{e)Y=^[^^j^ mx^dt-^{e)yFT. 

Next, we note that (see [2]) 



Hence, we have 

cT 



^j JEX?dt-n0) = ^l (]EX2-IEC/2)dt 

Thus, by Proposition 12.21 we obtain that as T tends to infinity 

(3.4) ^{tIo ^^'"^^ " "^^^0 ^ ° 

Therefore, it suffices to show that 

Ft ^AA(0,o-2), as T ^ oo. 
Now by Lemmas lA.H and lA.2l in the Appendix A, we have 

||L>,Ft||| -^ 2a^ and IE(F|) 
So the result follows by applying Proposition 12.11 
Now we are ready to state the main result of this section. 






D 
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Theorem 3.2. Assume we observe Xt at discrete points {t^ = kA^jk = 
0,1,..., A^} and T/v = A^A^v- Assume we have An — )• 0, T/v — ?• oo and 
NA'j^ -^ as N tends to infinity. Put 



1 ^ 
(3.5) /22,7V = 7^ V ^t^ Atfc and 0^ := * ^ (/Ja.Tv) 



Then 9 is a strongly consistent estimator of drift parameter 9 in the sense 
that as N tends to infinity, we have 



(3.6) 9n ^9 

almost surely, where "^^^ is the inverse of function ^ given in i2.13\) . More- 
over, as N tends to infinity, we have 



(3.7) ^TNi9N-9)^M{0,ai) 

where 



2 

and a^ is given by 113. 1\} . 



For the proof of Theorem l3.2l we need the fohowing version of strong law of 
large numbers. For the proof, see |13] . 



Lemma 3.2. Let Zj. be a sequence of random variables such that 

N N 

(3.9) ^j;]E[ZfeZ,]<CiV^ 

k=i j=i 

for some constant C < oo and /3 S (1, oo). Then for any 7 > f 

1 ^ 
fc=i 

almost surely as N tends to infinity. 
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Proof of Theorem \3.S[ Applying Lemma 12.41 we obtain 

rTpf 



13 



Ji2,N -J^ I ^t'dt 



1 



/Tn 



k=l 



E iK-^?)^t 



N 



^^^ E ^^p i^«i / \^t,-xt\dt 



/Tj 
2Yi{H 



< --'Z1-' I J2 sup \Xu\ j {tk- t)dt 



'Ti 



N 



^;^ue[tfe_i,ifc] 



+ 



2 



Y^ sup \Xu\Y2{H,e,[tk^i,tk]) r {tk-t)dt\ 

^_]^ MG[tfc_l,tfe] Jtk-1 J 

2^ sup |X„|y3(^,^,[ifc-i,ifc]) / (ife-i)dM 



N 



=: /i + /2 + /a + /4. 

We begin with I^. Define the random sequence {Zk}k>i as follows 

Zk= sup \Xu\Y4,{H,e,[tk-i,tk]). 
Me[tfe_i,tfc] 

Then using Cauchy-Schwartz inequality and Lemma 12.41 we obtain that for 
every j and k 



1 1 

4\ 4 /Tipr „,,„ I V ii4\ 4 



lE[ZfcZ,]< (1E[ sup |X„|]4jMlE[ sup |X„ 



<c 



where C is a constant independent of j and k. Moreover, we can take C to 
be independent of A^, because of the fact 

sup \Xu\ < sup \Xu\ 
ue[to,ti] ue[o,i] 

and 

Y4{H,e,[to,h])<Y4{H,e,[0,l]). 

Now we apply Lemma 13.21 with /3 = 2 to obtain that 

1 ^ 
— -V sup \Xu\Y^{H,e,[tk-i,tk])^0 
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almost surely for any 7 > 1. Take now 7 such that 1 < 7 < | H ^ with 

an enough small e > 0. With this choice, we have 

2 1 ^ 

1 ^ 
-y^A^-'N^-'—Y, sup \Xu\Y4H,e,[tk-i,tk]) 
i iV ' f — ' „cr+, , f, 1 



-^0 



almost surely, because the condition A^A^ — )• implies that 



N 

1 



/ 2H+l-2e \ 7-9 _l 

Treating Ji, /2, and /s in a similar way, we deduce that 



0. 



(3.11) y/TN 



M2,Ar - 7^ , 
J^N Jo 



1 r^^ 

x?dt 







almost surely. Moreover, we have convergence (|3.6|) by Lemma 12. 3[ To 
conclude the proof, we set fi = "^{9) and use Taylor's theorem to obtain 

TV (dj^ -e)= A^-i(^)y^(^2,jv - *W) 






for some reminder function Ri{x) such that Ri{x) — )• when x — )• ^{9)- 
Now, by continuity of — ^^^ and ^^^, we have that Ri is also continuous. 
Hence the result follows by using p. lip . Theorem 13. H Slutsky's theorem 
and the fact that 

T,""^^^ = W(9-y 

D 

Remark 3.1. Note that we obtain a consistent estimator which depends on 
the inverse of function ^. Despite we proved that such inverse exists, up to 
out knowledge there exists not an explicit formula for the inverse. Hence, 
the inverse must be computed numerically. 

Remark 3.2. The Theorem \3.2\ imposes different conditions on the mesh 
Atv to obtain strong consistency of the estimator 9. One possible choice for 
the mesh satisfying in such conditions is A^r = °^ . 

Remark 3.3. Note that we obtained strong consistency of the estimator 
9 without assuming uniform discretization of the partition. The uniform 
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discretization will play role in estimating the Hurst parameter H and for 
analysis of the estimator 6, see sections\^ and\5^ 

4. Estimation of Hurst parameter H 

There are different approaches to estimate Hurst parameter H of fractional 
processes. Among all, here we consider an approach which is based on fil- 
tering of discrete observations of a sample path of process. For more details 
we refer to [8j and [6]. 

Let a = (ao, ai, • • • , ai) G IR^^^ be a filter of length L + 1, L G IN, and of 
order p > 1, i.e. for all indices < q < p, 

L L 

Y,^jf = and Y^ajfy^O. 

j=0 j=0 

We define the dilated filter a^ associated to filter a by 

2 I Ofc', k = 2k' 
I 0, otherwise 

for < A; < 2L. Assume we observe Xt at discrete points {tk = kA]\f,k = 
1, . . . , A^}, such that the mesh Ajv — s- as A^ tends to infinity. Denote the 
generalized quadratic variations associated to filter a by 

(\ 2 
L 
j=0 

We consider the estimator Hjy given by 

(4.1) ^iv = i log/^''^' 



2 °^ Vn. 



a 



Assumption (A): 

We assume that for any real number r such that < r < 2p and r is not an 

even integer 

L L 
i=0 j=0 

Example 4.1. A typical example of a filter with finite order satisfying in 
assumption (A) is a= (1, —2, 1) with p = 2. 

Theorem 4.1. Let a be a filter of order p > 2 satisfying in assumption 
(A), and the mesh Ajv = N~°^ for some a G (^, jjj^)- Then 

H]y — > H 

almost surely as N tends to infinity. Moreover, as N tends to infinity, we 
have 

^/N{Hn - H)) A J\f{0, r{H, 6, a)) 
where the variance T depends on H , 6 and the filter a and is explicitly 
computed in [61 and also given in |4j. 
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Remark 4.1. It is worth to mention that when H < |, it is not necessary to 
assume that the observation window T^ = NA]\f tends to infinity, whereas 
when H > ^, we have to have T]\j- tends to infinity, see [8j. Notice that 
H >\ if and only if jjjz^ < 1 • 

Proof. Let vu denote the variogram of the process vjj . By Lemma 12. H we 
have 

vu = Ht^^ + r{t) 

as i — >• where r{t) = o{t ). Moreover, we have that r(t) is difFerentiable 
and direct calculations shows that 

Hence the claim follows by following the proof in ^J for the fractional 
Ornstein-Uhlenbeck process of the first kind and applying the results of Istas 
and Lang [8]. To conclude, we note that the given variance is computed in 
(iSj, p. 223). 

D 

5. Estimation of drift parameter when H is unknown 

Instead of ^{9), we consider ^{0,H) to take into account the dependence 
on H. Let ^ = ^{9,H) and define a new function 

Then implicit function theorem implies that there exists continuously differ- 
entiable function g{fj,, H) such that 

where 6 is the unique solution to /i = ^{9,H). So for every fixed H, we 
have 

Moreover, by chain rule we obtain 

^ = ^9{mH),H) = — + - — . 

Here -^ and qjj are known and hence we can compute qjj . Let H^ be given 
by ()4.ip for some filter a of order p > 2 satisfying assumption (A) and Jl2,N 
by ()3.5p . We consider the estimator 



(5.1) On = g{V'2,N,HM). 

Now with all above assumptions and the notations, we have the following 
result. 

Theorem 5.1. Assume that the mesh A^r = A^""^ for some number a G 
(2, 4ffT2 ^^)- Then the estimator 9^ given by (j5.ip is strongly consistent, 
i.e. as N tends to infinity, we have 

(5.2) 9n ^9 
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almost surely. Moreover the following central limit theorem 



(5.3) ^(eN-9)^N{<d,al) 

holds as N tends to infinity, where variance a^ is given by (37c 
Proof. First note that 

Tn [On ~^) = VTn i9{'fj'2,N,HN) - g{V'2,N, H) 
+ \/TN[9(fi2,N, H) - g{fj,, H) 
Now convergence 

y/f^{g{fi2,N,H)-g{ii,H)) AAA(0,ai) 
is in fact Theorem 13. 2[ Moreover, we get by Taylor's theorem that 



(5.4) 



Tn (g{li2,N , H]sf)-g{'fi2,N , H)) = -^{V'2,n,H)^JTn{Hn - H) 



+ -prT^{V'2,N , H)R2{'jl2,N ■, Hn)\/Tm{Hm — H) 

for some reminder function i?2 which converges to zero as {fi2,N,HN) — ^ 
{fi,H). Therefore, by continuity and Theorem 14.11 as N tends to infinity, 
we obtain 

V^N [g{V'2,N, Hn) - g(jj.2,N, H)j > 

in probabihty. Hence, using Slutsky's theorem we obtain 

AA(0,ai). 

by continuous 

D 



Tn {On 



To conclude the proof, we obtain ()5.2p from equation ([5 
mapping theorem. 
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Appendix A. Computations used in the paper 
Lemma A.l. For Ft given by i3. 2\) we have 

(A.l) 

Proof. It is sufficient to show that 

(A.2) IE \\D.,FtK-, -1E\\D.,FtK-, -^ 0. 



I n 7? I|2 JlL 0^2 



\DsFt\\}^-1E\\DsFt\\1 



This implies 



\D,Ft 






M\DsFt\ 



n 



21E{Fr^ 



T)^ 



which in turn implies (jA.ip . Now we have 



DsFt 



inaisr))- 
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Hence using the Remark 12.21 we can write 

t2H-2 f-T f-T 



DsFtWI = '"^f / / /fQ(n,.))/fQK-)) 



eH — eH 



2H-2 



dvdu. 



Using multiplicative formula for multiple Wiener integrals, we see that 

lH~g{u,-))I?{Kv,-)) 

= {giu, ■),g{v, ■))^ + /f {g{u, ■)g{v, •)) 

=:Ai + A2. 

Here Ai is deterministic and A2 has expectation zero. Hence, in order to 
have ()A.2p . we need to show that 

T Jo Jo 
Using Minkowski's inequality for integrals, we get 



(A.3) 



IE 



dvdu 



0. 



'E 



1 
T 

< 



Jo 

T rT 



eH — en 



2H-2 



dvdu 



^ Jo Jo 



eH — eH 



2H-2 



dvdu. 



Notice that A2 is a double Wiener integral with respect to process G. Hence 
using isometry property, we obtain 

W.{A2f = 2\\~g{u,rg{v,-)f 



2alH^''-^ 



g{u, ui)g{v, vi)g{u, U2)g{v, V2) 



[0,T]* 



^ ^{Ul+U2+Vl+V2)[jj-1) 



U2 lii 



2H-2 



V2 1^1 

e H" — e"H" 



2H-2 



duidfidn2du2 



where the function (/)(u) is defined by 

rT rT 



U 



JO 



5(^,zi)^(n,zi)e(^^+^^)(i-i) 



£2 £1 

eH — eH 



2H-2 



dzidz2- 



So, it remains to show that 



(A.4) 



T rT 



JO 



V</'(n)(^(w)e("+^Hi-i) 



eH — eH 



2H-2 



dvdu — )■ 0. 



We begin by showing that (f){u) is bounded as a function of u. We denote 
by C a constant which may vary from line to line. First notice that 



5(x,y) = ^' 



e -\^-y\ _^-d{2T-x-y) 



< c{e). 
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Taking account the symmetry and using change of variable z = e » , we 
obtain 



9 



(n) < C{e) / ~g{u,y) / e^^+^K^-i) [el - gf ) 



2H~2 



dxdy 



e'lT 



<C{0)B{l-H,2H-l) f ~g{u,y)dy. 

Jo 



It is straightforward to see that 

rT 



I 

Jo 



g{u,y)dy<C{e) 



where the constant is independent of T and u. Hence we have (j){u) < C{9) 
uniformly in u. Therefore, we obtain 

1 r'^ f-T , , . ., ^ -.^ „ „ 2H-2 



JO 



yJ<P{u)<l){v)e^''+<Ti-'^) 



en — ea 



eH — eH 



r'vT pu pT rwT pT pu 

/ + / / + / / 

/o Jo JVtJo JVtJVt 



dvdu 

2H-2 



~ T 
=: J1 + J2 + Js- 
For the term Ji , we obtain by change of variable that 



dvdu 
(«+-)(ir-i)(et-e^)'""'dz;dn 



2H-2 



J^ 



y/T /.I 



_ 1 

C 



z-^(l-z) 

e"7I 



2H-2 



dzdu 



as T tends to infinity. For the term J2, we obtain by Tonelli's theorem and 
change of variable that 



j, = l 






C 

Vt 



^-Hn _ ,fH-2^z^^ 



Similarly, for the term J3, we get 



T Jy/T Jy/T 



(-+^)iir-^)(el-elY''~'dvdu 



< C 



Vt 



T 



,-H(l _ zfH-2^^ ^ o_ 



Combine the estimates for the terms Ji, J2 and J3 together with the fact 
that the function (f){u) is bounded, we deduce (|A.4p . The proof is completed 
since we have implications (|A.4p =^ (|A.3p ^ (|A.2p =4> ()A.ip . D 
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Lemma A. 2. For Ft given by i3. 2\) . as T tends to infinity, we have 
(A.5) ]E[F|] -^ a^. 

Proof. First using isometry property, we obtain 



2It 



^[^t] - j;\\9\\'fi02 -■ -j7 



where 



It =alH^''-^ 



[0,T]4 



l{ui,viYg{u2, t,2)e(^-i)("^+''^+"2+^2) 



'"■2 "1 

e"H — e~H 



2H-2 



V2 n 

QH — eH 



2H-2 



dnidu2dwidv2. 



Recall that 



g{x,y) = —e 



_}_p-e\x-y\ _ J_ e{2T-x-y) 

29 26 



We first show that we can omit the second term ^e '^ ^ '^' in the function 
g. To see this, we have 

-e{2T-u^~vi)~(^^^^^^^^{jj^l){m+v-,+U2+V2) 



[0,T]' 



U2 ■"-[^ 



2H-2 



V2_ V2_ 

eH — eH 



2H-2 



dnidu2dfidti2 



< C{e) I e 

'[0,T]4 



-e{2T-u-i-vi) (jj-l){ux+vx+U2+V2) 



X 


^2 

eH - 


- eH 


2H-2 


= 


C{6) 


[f 

Jo 


Jo 



V2_ V2_ 

eH — eH 



2H-2 



dnidu2dwidu2 



e{T-v)+{jj-l){v+u) 



eH — eH 



2H-2 



dvdu 



By change of variables v = T — v,u = T — u, and then x = e h ^ y = e h 
this is same as 

[\ r,x(^-i)v^iy-xp^-2dW 

Je H Je H 

Let now x < y. By change of variables 2; = -, we obtain 
"1 ry 



[ , r,x(^-i)V^l2/-x|2^-2dxdy 

Je~TT Je'TT 



< 



Jo Jo 



<-^B{{e-l)H + l,2H-l) 

which converges to zero when divided with T and let T tends to infinity. 
The case x > y, can be treated in a similar way. Hence it is sufficient to 
consider the function 

J_p-(>\^-y\ 
26 
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instead of g{x, y). We shall use L'Hopital's rule to compute the limit. Taking 
derivative with respect to T, we obtain 






-e\T-ui\ -e\u2-V2\ {jj~l){T+ui+U2+V2) 



[0,T]3 



1^2 ^1 



2H~2 



2L "1 

eH — eH 



2H-2 



dnidii2di;idf2. 



By change of variables x = T — ui, y = T — U2 and z = T — wi , this reduces 
to 



diT ajrH^^-^ 



dT 



02 



-fe -e|j,-2| (l--^)(a;+y+^) 



X ( 1 — e H 
Therefore, we have 



[0,T]3 
2H-2 



e e 



e H — e H 



2H-2 



dzdxdy. 



lim ,^ 

T^oo dT 



diT _ alH^^~^ 



02 



X 1 - e H 



-ex^e\y~z\ {l~jj){x+y+z) 



e e 



[0,oo)3 
2H-2 



e H — e H 



2H-2 



dzdxdy. 



We end the proof by showing that the later triple integral, denoted by /, is 
finite. Use the obvious bound e~ '^~^' < 1, we infer that 



I < I e "-"e 

'[0,oo)3 

y \ Zll 2 

X ( 1 — e~H 



^^e(^~^)(^+^+'') 



e H — e H 

2H-2 



2H-2 



dzdxdy 



e{^-7i)y (l-e-H'\ dy 



oo roo 



ex {i^jj)(x+z) 



e e 



Jo 



e H — e H 



2H-2 



dzdx 



For the term /i , we obtain by change of variable u = e « that 

•1 



~H 



2//-2, 



Ii=C u~''{l-uy"-'du<oo. 
Jo 

For the term I2, we obtain by change of variables u = e~'H and v = e~Ti 
that 

h = C C C u^'-^^"v-''\u - H'^-'d^d^ 

JO JO 

+ / / ^(^-i)^^-^b-H'^-'dz;d« 



Jo Jo Ju 
= h,! + -^2,2- 

For the term /2,i; we obtain by change of variable z = - that 

hi = C I J^'^ I z'^{l - zf^-^dzdu = -^B{1 -H,2H- 1) 
Jo Jo (^H 
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Similarly for the term /2 2, we get by change of variable z = - that 



V 



[2 

[3; 
[4; 
[5; 

[6; 
ir 

[8 
[9 

[10: 
[11 

[12 

[13 

[14 

[15 
[16 
[17 
[18 
[19 



l2,2 = -^Bi{e-l)H+l,2H-l). 
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